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The conductance and the Fano factor in a graphene sheet in the ballistic regime are calculated. 
The electrostatic potential in the sheet is modeled by a trapezoid barrier, which allows to use the 
exact solution of the Dirac equation in a uniform electric field in the slope areas (the two lateral sides 
of the trapezoid). A special attention is devoted to asymmetry with respect to the sign of the gate 
voltage, which is connected with the difference between the Klein tunneling and the over-barrier 
reflection. The comparison of the developed theory with the experiment supports the conclusion 
that the Klein tunneling was revealed experimentally. 

PACS numbers: 73.23.Ad,73.50.Td,73.63.-b 

I. INTRODUCTION 

The Klein tunneling is one of the most important manifestations of the relativistic Dirac spectrum in graphene 2 . 
In this process an electron crosses a gap between two bands, which is a classically forbidden area, transforming from 
an electron to a hole, or vise versa. The Klein tunneling was well known in the theory of narrow-gap semiconductors 
under the name of interband or Landau-Zener tunneling. In the framework of the band theory the electron wave 
function for the states close to a narrow gap between two bro ad bands must satisfy the Dirac-like equation 3 -*. This can 
be demonstrated within the model of nearly free electrons^'. So the analogy with the relativistic electrodynamics was 
well known and exploited in the theory of semiconductors. For example, Aronov and Pikus^used the pseudo-Lorentz 
transformation (with the Fermi velocity playing the role of the light speed) treating the effect of the magnetic field 
on the interband (Klein-Landau-Zener) tunneling. This method was used by Shytov et a?!^ for graphene. 

One may expect to reveal evidence of the Klein tunneling from observations of char ge trans port and shot noise in a 
graphene sheet in the ballistic regime, which are now intensively studied experimentally^* 10 * 11 *. Analyzing conductance 
and shot noise in a ballistic graphene sheet a commonly accepted assumption was that under electrodes the graphene 
is strongly doped. A further assumption, which simplified a theoretical analysis, was that the level of doping changed 
abruptly. This led to a rectangular potential barrier for electrons in a graphene sheetP^LL^LLiJ. One may expect that 
in reality the doping level shoul d vary c ontinuously. Smooth finite-slope potential steps were analyzed theoretically 
for p — n transitions in graphen d 15 * 16 * 17 *. A possible model for the potential barrier might be a trapezoid s hown in 
Fig. [l] Recently experimental investigations of transport through tunable potential barriers were reporte d 18 * 19 * 20 !, 
which focused on observed asymmetry of the dependence of resistance on gate voltage with respect to the sign of 
voltage measured from the electrostatic potential of the Dirac point in the sheet center. In the diffusive regime 
asymmetry was attributed to scattering by charged impuritie d 21 * 22 *. In the ballistic regime asymmetry is related with 
the Klein tunneling 20 . 

Qualitatively the origin of asymmetry in the ballistic regime is illustrated in Fig. [T] We consider the limit of a 
small voltage bias, i.e., a difference of electrochemical potentials in leads is very small. Then only the states near the 
Fermi level contribute to the conductance. The small voltage bias drives electrons to the left and holes to the right. 
If the gate voltage V g is positive (Fig.[lji), the Fermi level crosses only the conductance band of electrons. There is no 
classically forbidden zone for electrons moving above the barrier, and the transmission of the sheet is restricted only by 
the over-barrier reflection. On the other hand, if the gate voltage V g is negative (Fig. [lb), the Fermi level crosses the 
conductance band in the electrodes and the valence band in the sheet. The sheet is classically forbidden for electrons, 
and the charge transport is realized via holes, which annihilate with electrons moving from the left to the area of 
the slope of the width d. This is the process of the Klein tunneling. In the limit of a very steep slope (rectangular 
barrier) the probability of over-barrier reflection exactly coincides with the reflection from the band b oundar jJ12U2l In 
general these probabilities are different: the probability of over-barrier reflection decreases with the growth of positive 
V g , while at large negative V g the Klein tunneling essentially restricts transmission and reflection probability remains 
finite. Therefore experimental detection of asymmetry provides evidence of the Klein tunneling^. 

For studying transport through the trapezoid barrier one should find the electron wave function inside the slope 
areas of the width d, where the electron is subject to a uniform electric field. The Dirac equation in a uniform 
electric field has an exact solution in terms of confluent hypergeometric functions, which was found long time ago by 
Sauter^. He used it calculating the probability of the Klein tunneling. Kane and BlountP gave an exact solution of 
the Dirac equation in terms of Weber parabolic cylinder functions , which are directly connected with the confluent 
hypergeometric functions^ used by Sauter 23 . Cheianov and Falkc^ also found Sauter's solution but by the method 
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FIG. 1: (color online) Trapezoidal electrostatic barrier in a graphene sheet. The thick broken solid line shows the electrostatic 
potential for the case, when the Fermi level in the sheet center crosses the Dirac point, (a) A positive gate voltage V g shifts 
the Fermi level up to the conductance band. The electron at the Fermi level is crossing the sheet from left to right without 
leaving the conductance band, (b) A negative gate voltage V g shifts the Fermi level down to the valence band. The electron 
crossing the sheet must tunnel from the conductance to the valence band on the left (inside the circle) and tunnel back to the 
conductance band on the right. Reflected waves are not shown. 



valid only for not very steep slopes while the solution is exact for any slope. Sauter's solution was also used for 
studying p - n junctions in carbon nanotubes . This solution will be an essential component of the present analysis 
of conductance and short noise in graphene, though exact solutions are known also for other types potential barriers, 
e.g. a barrier with exponential variation of the electrostatic potential. The paper presents calculations of the 
conductance and the Fano factor for a trapezoid potential barrier in a ballistic graphene sheet as functions of the gate 
voltage. 

Section [IT] presents the Dirac equation for graphene and the semiclassical analysis of the Klein tunneling. Section III 
analyzes the Klein tunneling using the known exact solution of the Dirac equation in a uniform electric field. Section IV 



studies transmission and reflection of electrons propagating across the trapezoid barrier formed in the graphene sheet 
by doping. Its results are used in Sec. [Vjfor calculation of the conductance and the Fano factor of the sheet. The last 
Sec. |VI|is devoted to comparison with experiment and concluding discussion. 
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II. DIRAC EQUATION AND SEMICLASSICAL ANALYSIS 

The Hamiltonian of the graphene in the presence of the electrostatic potential V(x), which depends only on the 
coordinate x, is 

H = v F (a x p x + a y p y ) + eV(x)I, (1) 

where vp is the Fermi velocity, p x = —iTid/dx and p y = —ihd/dy are components of the momentum operator, / is a 
unit 2x2 matrix, and a x , a y , and a z are Pauli matrices of the pseudospin. The eigenstates are spinors 

•<*•»>-($)■ < 2 > 

where the components ipn are amplitudes corresponding to the eigenvalues ±1/2 of the spin matrix a z . The compo- 
nents of the eigenstates with the energy e satisfy the equations 

~ ( + IT I = ^(^h 



dx dy 
.d_ _ d_ 
dx dy 



ip r = K(x)i/)i, (3) 



where K(x) = [e— eV{x)]/TiVF- For the further analysis of the exact solution it is more convenient to perform rotation 
by 90° around the axis y in the spin space introducing the spinor 



where ij)± — (±?/>j + ipi)/V% are amplitudes corresponding to the eigenstates of the spin matrix a x . After rotation 
the Hamiltonian becomes^ 

H = v F (a z p x + OyPy) + eV(x)I, (5) 
and the amplitudes ip± satisfy the equations 

- l ^~^y =]C{x) ^ 
.di>- dtfj+ 

,_ + _ = £(*),&_. (6) 

Before analyzing the exact solution of these equation for linear function )C(x) (see the next section) it is useful to 
present a less accurate but physically more transparent semiclassical analysis. The semiclassical solution of Eq. ^ 
for the x dependent potential is 

*(*.y) = ($!) = ( 4+^5- ) \fl/ rHx ' )dx ' +ik » v , (7) 
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where k(x) = \K(x)\ and k x (x) — ^JIC(x) 2 — ky. The x component of the current in our presentation is 

j x = ev F ^a^ = ev F {\^ + \ 2 - |^_| 2 ), (8) 

and the spinor is normalized to the current equal to j x — evp. 

One can use the semiclassical solution for the analysis of the intcrband transition, which takes place for the case of 
nee' ative V g (Fig. []J>). Fi gure pi shows an electron with the wave vector \n(k x ,k y ), which moves from the left to the 
right in a uniform electric field (/C = —ax). This is a zoom in on the Klein-Landau-Zener process inside the circle 
shown in the left-slope area in Fig. [T(d and is a slightly revised version of Fig. 114 in the book by ZimarP. At some 
moment the electron enters a classically forbidden area where the x component of the wave vector becomes imaginary: 

k x (x) — ±i^/fc^ — IC(x) 2 . The point where k(x) — is the turning point of the classical trajectory. In the point where 

K{x) = the electron crosses the border between states with positive and negative band energies (particle and hole 
bands). Let us call it "crossing point". Choosing the crossing point to be at x — 0, the classically forbidden area 
extends from x — —x c to x c — k y /a. Then according to the semiclassical approach the probability of the tunneling is 



(9) 

This semiclassical analysis is not expected to provide an accurate pre-exponcntial factor. But remarkably the estima- 
tion fully coincides with the exact result given below. 



A" 



exp 



k 2 — JC(x) 2 dx 



III. EXACT SOLUTION FOR AN UNIFORM ELECTRIC FIELD 

Assuming JC{x) = —a(x — xq) the general exact solution of Eq. ([6l is^ 



*(x v) - ( ^i-rK."-; -r^ains./V ) e ik v v (in) 



C 1 F(e, K ) + C 2 G(C,K) 
CiG(e,K) + C 2J F(e,K) 

where £ = (x — Xo)y/a, k = k y /y/a, xq is the coordinate of the crossing point where JC(x) = —a(x — xq) = 0, 

1 Ki 1 o\ , . /O t r f t Ztv 3 



F&k) = e"« 2 / 2 M (-^ , = -<e"^/ 2 M (l - ^, |,^ 2 ) , (11) 



and 



M{a,b,z) = J2^r Z J^ (a)o = l, (a)„ = a(a + l)(a + 2).„(a + n - 1) (12) 

71=0 

is the Kummer confluent hypergeometric function^ satisfying Rummer's equation 



The two constants C\ and Ci are determined by the boundary conditions. According to the known asymptotics of 
the Kummer functiomPH, at large distances £ — > ±oo 

k) = F 00 ( K )e-« 2 / 2 r 3 /2 ! G{t k) = _|_ Goo(K)e ^ 2 /2 rl « 2 /2 ; (14) 

where 

^oo(^- r(1/2 + m2/4) , Goo^j- 2r(i-«V4) ' ( j 

The functions F'oo(k) and Goo(k) satisfy the relational! 

|Foo(«)| 2 = - |Goo(K)| a = . (16) 
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Thus the exact solution at large distances reduces to a semiclassical solution since 

e i£ 2 /2£-iK 2 /2 ^ e « /J" k x (x')dx'+ik y y _ g J* a 1 x' 2 -k^dx' +ik v y 

One can see that very far from the crossing point JC(x) = (£ = 0) the electron moves parallel to the x axis, and the 
asymptotic of the exact solution at £ — > ±oo is 

* = ( ) e-^^^CCifoo T C 2 G^) T ( J) e^+^CdGoo T C 2 F* ). (18) 

We investigate the process shown in Fig. [2] an electron moving from left either transforms after the Klein tunneling 
to the electron with negative energy (a hole moving from the right to the left), or is reflected backward. Thus at 
x — > oo (£ 3> re), the solution should transform to the semiclassical solution 

»(*,») - ( i^+S, ) ^-e-'i:^'^'^ _> (J ) e-^/*^". (19) 

Using Eq. (18 1 and the identity l^ool 2 - l G oo| 2 = 1 

one obtains 

Gi=i^, G 2 = G oc . (20) 



Then the asymptotic at £ — ► — oo is 

1 fl 



( J ) e -« a /a+iM + ^ ( J ) e'^/a-Hfc,^ (21) 



where tff and are amplitudes of transmission and reflection determined from relations 

^- = \F 00 \ 2 + \G 00 \ 2 ^e~^ 2 / 2 , r JL = 2 F* 00 G 00 . (22) 

This yields the exact probability of the Klein tunneling Tk = \tx\ 2 = eT 7 "^ ' , which coincides with the semiclassical 
result Eq. 



IV. TRANSPORT ACROSS THE TRAPEZOID BARRIER 



A. Scattering at the left side of the barrier 

Now let us consider electrons moving across the trapezoid barrier shown in Fig. [l] In the area of the slope [x < 0) 
electrons arc in a uniform electric field, while the area x > is field-free. We look for a solution, which in the field-free 
area x > is a plane wave with the current evp. For positive gate voltage V g the plane wave corresponds to an 
electron of positive energy [JC(x) > 0] 

Ci k x +ik v \ rrr 
V , kJUk, J ^J-e^+^y, (23) 

while for negative V g the electron has a negative energy, and its group velocity x component has a direction opposite 
to the x component of the wave vector k(/c x , k y ): 

V) = ( A +tsL ) ^Ee-^+^vV. (24) 

The constants C\ and G 2 in the exact solution for x < are now determined from the continuity of the spinor at 
x = 0. Introducing the reduced gate voltage v = eV g /TiVF\/a, the fitting point x — corresponds to the argument 
£ = — v of the Sauter's solution, which is negative for V g > (over-barrier reflection) and positive for V g < (Klein 
tunneling). From fitting one obtains the following expressions for C\ and G 2 valid for the both signs of V g : 

k + k x + isignw k v . k-k x - isignu k y 

1 = 2k ^' K ' 2k ^' ^' 

fe - k x - isignw fcj, fc + k x + isignu fc y 
G 2 = — y -F(v,K) — y -G(v,n). (25) 
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Here the properties F(—£, re) = F(£, re) and G(— £, ft) = — G(£, re) were used. 

Using Eq. ( 18 1 and the calculated values of C\ and G2 one obtains the following asymptotic expression at x 



*=i(j) e -<«»+«.. + ri(:)e' 



C/2+ik y y 



where the amplitudes of the transmission, t\, and of the reflection, r 1; are given by 

1 k + k x + isignu k v . k — k x — isignw k v 

— = — — " nil) KA — 

t>\ 2^/ kk x 



v K) _ l ^""" Giv, k), 



2~\/ kk^ 



-co: 



(26) 



(27) 



ti 



isignv k y 



isignw k y 



2 \f k \\jrr- 



F(v, re) 



Here 



F(v, re) = F*(v, K)i?oo(«) + G{v, k)Goo(k)*, G(v, re) = G>, re)^ + reJG^re)* 
Eventually the transmission probability at the left slope of the trapezoid barrier is given by 



k ~\~ ki 



\G(v, k)\ 2 - signw-^Im{J r (u, n)*Q(v, k)} 



(28) 



(29) 



(30) 



2 kj^ 2 kj^ 

According to Fig. [l]the electric field is absent deep in the electrode at x < — d . So in the point x = — d the field 
solution must transform to a plane-wave solution again. However, there is no significant reflection of the electron 
in this area, if the asymptotic expression Eq. (26 1 is valid at x = —d. Indeed, according to Eq. (26 1 the electron 
propagates normally to the barrier and cannot be reflected. The conditions for it are d)| ~ dy/a ^ 1 and 
d)\ k. In particular, this means that 

(31) 

where fco = eVo/hvpy/a is the modulus of the wave vector at the Fermi level inside the electrode at x < —d. As far as 
this condition is satisfied, neither k$ nor d affect the results of the analysis, which depend only on their ratio a = ko/d 
proportional to the electric field in the area of the slope — d < x < 0. 

Let us consider various limits of the obtained expression for the transmission. The very steep slope (rectangular 
barrier) corresponds to very small \v\ and n. In this limit F(v,k) = F*(v,k) — ► 1 and Q{v,n) = G*(v,k) — > 0, and 
one obtains that transmission independent of the sign of V g . i.e., the difference between the Klein tunneling and the 
over-barrier reflection vanishes: 



Ti 



k + k x 



(32) 



This result is valid for small \V g \ much less than TiVFy/a/e (\v\ <C 1). 
Let us consider now the opposite limit of high \V g \ when \v\ 3> re, and 

F{v,k) -» [|F 00 (re)| 2 -sign t ;|G 00 (re)| 2 ]e™ 2 / 2 , Q{v, re) -» 20(-v)F oo (re)G 00 (re)* e - i! ' 2 / 2 . 
According to the definitions v and re the condition \v\ 3> re means that k w k x ^> k v , i.e., the incident electron moves 



(33) 



nearly normally to the barrier. Then one should keep only the first term in Eq. (30 1 

1 4 



Ti = 



[|^ co (re)| 2 -sign«|G co (re)| 2 



[e™ 2 /2 + l - sign-y ( e ™ 2 /2 _ 



(34) 



This yields the probability of the Klein tunneling given by Eq. ([9 ) for negative v and ideal transmission T\ = 1 for 
positive v, i.e., for large positive v the over-barrier reflection vanishes, 

For further calculation of the transmission of the whole barrier one needs also to know the parameters of the process 
time-reversed with respect to scattering at the left slope of the barrier. The time-reversed state corresponds to the 
negative current —evp (from the right to the left) and is described in the field- free region x > by the spinors in 
Eqs. (23 1 and (24) after replacing k x by —k x . Also the roles of an incident and a reflected wave in the asymptotic 



expression Eq. (26 1 are interchanged, and the transmission and the reflection amplitudes are determined from those 



for the original process as 



ti (k Xl ky) 



r i{ k x ,ky) fi{k x ,ky) 

tl(~k x ,ky) tl(k X ,ky) 



ti( k x , ky) 



(35) 
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B. Scattering at the right slope of the barrier 

The analysis of this process is similar to that done for scattering at the left slope: one should fit the exact solution 
in the region x > L of the uniform electric field, which has now an opposite direction with respect to that at the left 
slope, i.e., K.(x) = a(x — xq). Without repeating all details we summarize here the results. 

The exact solution at the right slope is complex conjugate to that at the left slope. Asymptotically the exact 
solution in the field region is a wave propagating to x — > oo: 



(36) 



The solution in the field area x > L must fit to the solution in the field-free region x < L, where there is a superposition 
of an incident and a reflected wave, which is 



for positive V g (over-barrier reflection) and 



h, r ( 1 I •%X~\~iky 

JKx+ik y y + I 2 + k 

*2 I ~2 + k 



-ik x x-\-ik y y , 

^2 



1 _i k x iky 

2 1 X~" 

1 _|_ fox IK 

~2 k 



for negative V g (Klein tunneling). The fitting at the point x = L gives 



1 

h 
t 2 



k + k x — isignu k. 



2 "\J kkfjQ 

k + k x + isignv k y 



k — k x + isignu fc. 



Q(v,k) + 



2 k}\t<£ 

k — k x + isignu fc. 



^(w,k) 



^ik x x-\-iky y 



T(v, k) 



isignu fe x L 



D isign« 



2\/ kk% t ^'\j J\iJ\ij^ 

This yields the transmission probability T2 = 1^2 1 2 equal to T\ at the left slope of the barrier. 



(37) 



(38) 



(39) 



C. Transmission through the whole barrier 



The transmission through the whole barrier depends on whether the electron can pro pagate inside the barrier 
coherently without losing its original phase^ [see the phase factors e ±tS1 g nvk * L m Eq. ( 39 1] . This is not the case for 
L long enough. Then the total transmission T is combined not from amplitudes but from probabilities, treating two 
slopes as uncorrelated scatters. Keeping in mind that T\ — T2 one obtains^ 



T = 



TjT 2 
1 — R1R2 



2-Ti' 



(40) 



where = |^i,2| 2 = 1 — ^1,2 are probabilities of reflection at the left and the right slope. 

If the phase correlation takes place one should look for the solution in the electric field at x < 0, which at x 
has the same form as Eq. ( 26 1 , 



-( 1 

t \ o 



-if; /2+ik y y 



t 



jf/2+ik y y 



(41) 



but the transmission and the reflection amplitudes are determined now from fitting to the spinors Eqs. (37) or (38 1 
in the point x = 0. This yields the following expressions for them [compare with Eq. (9) in Ref. I16j: 



1 



1 



t\ti t\t<z 



cos k x L — — sin k x L 

fox 



J-(v, k) 



cos k x L 

2k, 



ik 
k x 



sin k. r L 



y 



sin k x LJ r (v : k)Q(vj k). 



(42) 



These expressions can also be used for the case, when there are no propagating modes in the field-free area < x < L, 
i.e., k x is imaginary and corresponds to an evanescent state in the classically forbidden barrier area. In this case 
one should analytically continue the expressions replacing smk x L/ik x with smhpL/p, where p = ik x is real if k x is 
imaginary. 
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FIG. 3: (color online) The plots (solid lines) reduced conductance g/go vs. gate voltage (a) and Fano factor F vs. gate voltage 
(b) for incoherent transport. The dashed lines show the results of simplified calculations approximating the transmission 
probability at negative V g with the probability of the Klein tunneling in a uniform electric field. The small black circle in (b) 
shows the value 0.125 of the Fano factor for vertical slopes (rectangular barrier) at voltages V g S> Tivp/eL. 



V. CONDUCTANCE AND SHOT NOISE 



A. Incoherent ballistic transport 



Let us consider the case of the graphene-sheet length of the L long enough, when the contribution of evanescent 
modes is not important, and there is no phase correlation between two tunneling events at the two sides of the barrier. 
Then the conductance does not depend on L and is given by 



9 



9o 



e\V g \/hv F 



Tdki, 



9o 



Tdn. 



(43) 



where go = Ae 2 W 
help of Eqs 



that W exceeds a 



a/irh, W is the width of the graphene sheet, and the transmission probability T is determined with 
(|30|) and ( 40 1 . One may replace summation over transversal components k y by integration assuming 



1 other spatial scales(L and 1/y/a). 
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Figure [3^, shows the reduced conductance 3/50 as a function of the reduced gate voltage v = eV g j \fahvp. At 
positive V g the conductance grows roughly linearly, which is related with the linear growth of the density of the 
states with the voltage. One may compare the conductance in Fig.[3]i with the conductance for a steep potential step 
(rectangular barrier): 

9r =90 I Tdk y =-g ^. ^ 







Here T = k x /k. The plot in Fig. |3p coincides with this dependence at very small v. Figure |3p shows the dependence 
of the Fano factor, 

fj*' T(l — T)d,K , 4 

F=- — n — - 45 

Jo Tdn 

on the reduced gate voltage v = eVJ yfahvp. The small black circle shows the value 0.125 of the Fano factor for 
vertical slopes (rectangular barrier)^ at voltages V g ^> fivp/eL. 

The left-hand parts of the plots (negative V g ) can be compared with the calculations assuming that the transmission 
T\ of a potential step is fully determined by the Klein-tunneling probability Tk [see Eq. (J9|], and T — Tk/(2 — T%) 
at any voltage. Then the conductance and the Fano factor are 

. . r\v\ 2{e^-lWK 

9 = 90 J Q F= rH *« • (46) 

JO 2e™ 2 -l 

They are shown in Figs. [3^, and [3Jd by dashed lines. This approximation is similar to that used by Cheianov and 
FalkoSl for a single p - n transition, but our analysis addresses two p - n transitions in series. One can see that at 
large negative gate voltage V g the conductance and the Fano factor are well described by the process of two sequential 
uncorrclatcd Klein tunnclings and saturate at the plateaus determined by 

g K = 0.4C% , F K = 0.329. (47) 

The subscript K stresses that these values are determined by the Klein tunneling only, and their observation is a 
direct evidence of the Klein tunneling. 

As mentioned above the content of this subsection refers to very large L. However this approach at some L large 
enough can fail because of scattering by disorder. Thus the validness of the approach is restricted with the window 
hvpjeV g <L<I where I is the mean- free path determined by disorder. 



B. Coherent ballistic transport 



In the coherent transport the total transmission amplitude and probability do depend on the length L of the field- 
free region. For finite length L the evanescent states participate in the transport, and in the expressions for the 



conductance and the Fano factor, Eqs. (43 1 and (45), one must repl ace the upper limit \v\ in the integrals by 00. 
Transmission T = \t\ 2 in these expressions is now determined by Eq. (42 1. The numerically calculated conductance 



and Fano factor as functions of v are shown for Ly/a = and L^Ja — 1 in Fig. [4] In the case of Lyfa — the barrier 
transforms from trapezoid to triangular. As well as in Fig. [3] dashed lines show dependences for the double Klein 
tunneling. One can see that at negative V g the conductance and the Fano factor for the coherent transport oscillate 
around the "Klein" plateau. 

The plots given above used the scales connected with the the length scale l/y/a, which is determined by the slope 
only. This length is more useful until it less than L, i.e, for not too steep slopes. However, in the opposite case of 
steep slope (a barrier close to rectangular) the length scale L might be more useful. Then it would be convenient to 
use Go = I/-R0 = 4:e 2 W/irhL (minimal conductance for a rectangular barrier) as a conductance scale and Tlvf/gL as 
a voltage scale. Figure [5] shows the plots of the reduced resistance R/Rq as a function of the dimensionless voltage 
vi = vL^/a = eVgL/fivp for three values of L^/a = 00, 10, and 5. Figure [6] shows the dependence of the odd part of 
the resistance -R oc y = [R(—v) — R(v))/2 on the reduced voltage Vi for two values of Ly/a = 10 and 5 (at Ly/a = 00 
the dependence is symmetric and -R oc y vanishes). At large voltages 2i? oc y oscillates around the Klein resistance 

= J_ = ^ h = (AS) 

K g K 1.612e 2 WV^ 1.612e 2 H/( 7 rn ) 1 /4 ' 1 > 

where no = k^/ir = eVo/hvp is the charger density in electrodes. 
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FIG. 4: (color online) The plots reduced conductance g/go vs. gate voltage (a) and Fano factor F vs. gate voltage (b) for 
coherent transport. Solid lines: 1 - L^/a = 0, 2 - L^/a — 1. The dashed lines show simplified calculations based on the 
probability of the Klein tunneling in a uniform electric field. 

VI. DISCUSSION AND COMPARISON WITH EXPERIMENT 

The odd part of the resistance was examined experimentally in Refs. 18 19 20 (see Figs. 3 in Refs. 18120 and Fig. 2 
in Ref. I19|l . There is some similarity between experimental curves and theoretical curves in Fig. [6] the majority 
of experimental curves also have plateaus around which the curves oscillate. This supports the claim of Stander et 
al. that they have found evidence of the Klein tunneling. The oscillations in the experiment look smaller and more 
chaotic than in the theory for coherent tunneling. This might be an effect of disorder. It is worthwhile to note that the 
appearance of Klein plateaus on theoretical curves is sensitive to the choice of the model: a plateau can appear under 
the assumption that the electrical field inside the transient area of the slope does not vary and the gate voltage V g 
does not approach too close to the voltage Vq, which determines the potential step (Fig. [TJ . The presence of plateaus 
on experimental curves demonstrates that this assumption is not so bad. For further comparison we use the data for 
the sample shown in Fig. 3a of Ref. [201 

Let us consider the dependence of the Klein resistance Rk (resistance at the plateau) on the electric field. If the 
width d of the transient area of the slope is fixed, the electric field is proportional to Vq, the latter being related 
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R/R 




FIG. 5: (color online) The plot resistance vs. gate voltage for coherent transport. The lines 1, 2, and 3 correspond to the values 
L^/a = oo, 10, and 5 respectively. 




FIG. 6: (color online) The plot odd resistance vs. 
values L^/a =10 and 5 respectively. 



gate voltage for coherent transport. The lines 1 and 2 correspond to the 



with the charge density no in the electrodes. In the experiment uq changes from 1.2 x 10 12 cm 



According to Eq. (48 1 Rk oc n ^ 4 must decrease in 1.4 times while in the experiment the plateau resistance decreases 
in about 1.6 times. A quantitative comparison of absolute values of the resistance is not straightforward because of 
the lack of information on the experimental value of the width d of the slope area (the width of the p - n transition) . 
Choosing the distance between the top gate and the graphene sheet as a rough estimation for d (34 nm for the sample 
under consideration), for no — 4.7 x 10 12 cm~ 2 Eq. (48) yields resistance 0.110 kf2 against about 0.125 kf2 in the 
experiment. 

To conclude, the paper presents calculations of the conductance and the Fano factor in a graphene sheet in the 
ballistic regime. The electrostatic potential in the sheet is modeled by a trapezoid barrier, which allows to use the 
exact solution of the Dirac equation in a uniform electric field in the slope areas (the two lateral sides of the trapezoid) . 
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A special attention is devoted to asymmetry with respect to the sign of the gate voltage, which is connected with the 
difference between the Klein tunneling and the over-barrier reflection. The asymptotic resistance for large negative gate 
voltage, when an electron crosses two p - n transitions in scries, is determined by the process of the Klein tunneling. 
The comparison of this asymptotic resistance (Klein resistance) with the experiment supports the conclusion that the 
Klein tunneling was revealed experimentally. 
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